home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / trap.f < prev    next >
Text File  |  2002-06-11  |  401b  |  20 lines

  1. C
  2. C PAGE 161: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: TRAP.FOR
  5. C
  6. C TRAPEZOID RULE PROGRAMMING EXPERIMENT FOR AN INTEGRAL
  7. C
  8.       DATA  N/60/ 
  9.       F(X) = 1.0/EXP(X*X) 
  10.       H = 1.0/(N-1) 
  11.       SUM = 0.5*(F(0.0) + F(1.0))     
  12.       DO 2 I=2,N-1
  13.         XI = (I-1)*H
  14.         SUM = SUM + F(XI)   
  15.    2  CONTINUE
  16.       SUM = H*SUM 
  17.       PRINT *,SUM 
  18.       STOP
  19.       END 
  20.